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We study Bragg spectroscopy of ultra-cold atoms in one-dimensional optical lattices as a method for probing 
the excitation spectrum in the Mott insulator phase, in particular the one particle-hole excitation band. Within 
the framework of perturbation theory we obtain an analytical expression for the dynamic structure factor S{q, uj) 
and use it to calculate the imparted energy which has shown to be a relevant observable in recent experiments. 
We test the accuracy of our approximations by comparing them with numerically exact solutions of the Bose- 
Hubbard model in restricted cases and establish the limits of validity of our linear response analysis. Finally 
we show that when the system is deep in the Mott insulator regime, its response to the Bragg perturbation is 
temperature dependent. We suggest that this dependence might be used as a tool to probe temperatures of order 
of the Mott gap. 

PACS numbers: 03.75.Hh,03.75.Lm 
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I. INTRODUCTION 

Recently, there has been a lot of experimental progress 
studying cold atoms confined in optical lattices. The defect 
free nature of the lattice potential, the long coherence times 
of the constituent atoms, and the high experimental control of 
the lattice parameters Qllj] make this a unique system for pre- 
cisely studying many-body physics. In particular, the exper- 
imental observation of the superfluid to Mott insulator quan- 
tum phase transition yl has stimulated much interest in this 
area of research. 

Perhaps one of the most important potential applications of 
the Mott insulator transition is to use it as a means to initial- 
ize a quantum computer registerr4l|5l|a,|3]. Particulary in this 
case, it is important to have tools for thoroughly character- 
izing the experimentally obtained Mott insulator states. The 
usual procedure for entering the Mott insulator regime is to 
begin with a magnetically trapped EEC (with almost all the 
atoms in the condensate), slowly load it into an optical lat- 
tice by increasing the lattice depth. One key piece of evidence 
for the quantum phase transition is the loss of global phase 
coherence of the matter wavefunction when the lattice depth 
increases beyond a critical value |3]. However, the loss of 
coherence could arise from many sources, such as the deco- 
herence induced by quantum or thermal depletion of the con- 
densate during the loading process |8] and therefore it is not 
a sufficient signature that the system is in the Mott insulator 
state. For this reason, in the experiments by Greiner et al.\y\, 
complementary evidence for the Mott insulator transition was 
provided by applying a potential gradient to the lattice to show 
the presence of a gap in the excitation spectrum. In this pa- 
per we show that Bragg spectroscopy, done by applying addi- 
tional laser beams independent from the lattice beams, is an 
experimental technique with the potential to thoroughly char- 
acterize the Mott phase. In addition to determining the energy 
gap, we show that it provides detailed information about the 
excitation spectrum, information unavailable using other tech- 
niques. Moreover, in contrast to applying a potential across 
the lattice Bragg spectroscopy is not susceptible to effects 



like Bloch oscillations and Zener tunneling. Furthermore, we 
show that different from the superfluid regime, the system's 
response to Bragg perturbation in the Mott regime is sensitive 
to finite temperature. This property might be used as a tool to 
probe temperatures of order of the Mott gap. 

Our analysis is based on a perturbative treatment, which we 
show to be applicable in the strong Mott regime that has been 
reached in current experiments [^. Although our approach is 
applicable only in the range of validity of first order pertur- 
bation theory, it has the advantage of properly including one- 
particle-hole correlations. Such correlations have dominant 
influence on the spectrum of the system and are not accounted 
for in mean-field treatments 1 10|. 

The organization of this paper is the following. In Sec. |ll| 
we introduce the basic formalism that describes Bragg spec- 
troscopy in an optical lattice and use a linear response ap- 
proach to calculate the imparted energy to the system. In 
Sec. |III]we derive the zero temperature response to Bragg 
spectroscopy of a translationally invariant lattice deep in the 
Mott regime and in Sec. lIVI we discuss the conditions required 
for our linear response analysis to be valid. In Sec.|y]we ex- 
tend the zero temperature analysis to finite temperature and 
finally in Sec. I VII we conclude. 



II. FORMALISM 

The typical Bragg spectroscopic procedure is to gently scat- 
ter atoms in an ultra-cold gas system with a moving potential 
of the form Vq cos{qx — uot). This type of experiment was 
first demonstrated by the NIST and MIT groups 1 11. .121 . In 
contrast to earlier experiments that used momentum as the re- 
sponse observable, here we choose to examine the imparted 
energy. In trapped systems this allows long excitation dura- 
tion that facilitates more precise spectral resolution. Energy- 
spectroscopy is not as well developed as momentum spec- 
troscopy but recent experiments have demonstrated the use of 
this technique lUsll . The Bragg potential is formed by the ac- 
Stark shift arising from a pair of interfering light fields (e.g. 



see fl41). In this paper we will always assume that the Bragg 
potential is generated independently of and is much weaker 
than the lattice potential. We therefore treat the scattering pro- 
cess with linear response theory. Using an independent set of 
beams to generate the Bragg potential also provides consid- 
erable flexibility in the range of q and lo values that can be 
obtained. 

In this work we consider one-dimensional bosons loaded 
in an optical lattice. Effective ID systems have been realized 
in recent experiments (for example refs. |9, 13]) by loading 
a Bose-Einstein condensate into a three-dimensional optical 
lattice, which is very deep in two directions. The dynam- 
ics is then restricted to the third, or axial, direction only. In 
this work we study the response of the system to Bragg per- 
turbation in the axial direction, assuming that the dynamics 
in the transverse directions is frozen. We consider a one- 
dimensional optical lattice which is sufficiently deep that the 
tight binding approximation is valid and assume that we can 
restrict the dynamics of the atoms to the lowest vibrational 
band. This applies when changing the lattice potential does 
not induce band excitations. This condition is fulfilled when 
the frequency, lo, of the Bragg perturbation is less than the gap 
between the first and second bands, and when the momentum 
transfer q is contained within the first Brillioun zone. A de- 
tailed analysis of the validity of this first-band approximation 
to study Bragg scattering of a dilute weakly-interacting gas 
in an optical lattice is found in 1 15] where the authors used a 
mean field approach combined with Bogoliubov analysis. In 
the single-band approximation and in absence of external po- 
tentials, the system is described by the Bose-Hubbard Hamil- 
tonian (BHH) ^ 
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Here a„ is the annihilation operator at site n which obeys the 
canonical commutation relations for bosons, J is the hopping 
matrix element between nearest neighbors, and U is the on- 
site repulsion energy. The sum {n, I) is taken over nearest 
neighbors. We use N for the total number of atoms and M for 
the total number of wells. In the tight binding approximation 
the Hamiltonian describing the Bragg perturbation reads 
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density 
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is defined 






operator 

y^M-l jn-m-'f ' iqmd where 

/ dxe^'^^(j)Q{x)(j)o{x — dn) is a geometrical factor that involves 
integration over Wannier functions (j)n(x), and d is the lattice 
constant. For deep lattices /" ex 5„ o lllSll . 

To analyze the Bragg spectrum of the system we study 
the energy transfer, which can be measured by time-of-flight 
techniques] 13]. Under linear response theory, the energy 
transfer is related to the so called dynamic structure factor, 
S{q, uj), which is given by 



S{q,uj) 
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where /q(w,y) = \{i\ Pq \j)\ 5(w-w,y), \i) and Ei are eigen- 
states and eigenenergies of the unperturbed Hamiltonian {I), 
g-PEi jg jjjg Qsuaj Boltzmann factor with (3 — l/ksT where 
ks is Boltzmann's constant and T the temperature, Z is the 
canonical partition function, and hujij = Ej — Ei. Because of 
the factor /^(cjy ), the system's response shows peaks when- 
ever the frequency of the Bragg perturbation matches the en- 
ergy difference between two eigenstates of the BHH. The peak 
height is proportional to the the transition probability between 
the two eigenstates, \{i\ pq |j)| . 

The total energy transfer after applying the Bragg perturba- 
tion can be shown to be given by 
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dt 
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where T.p is the duration of the perturbation and x{Q^ ^) = 
S{q,LLj) — S{—q, —uj)- Here we derive analytic expressions 
for the dynamic structure factor assuming we are deep in the 
Mott insulator regime, where treating the hopping term in the 
Hamiltonian as a perturbation is justified. 



III. ZERO TEMPERATURE RESPONSE 

In this work we assume a commensurately filled lattice with 
no external confinement, filling factor N/M — g and periodic 
boundary conditions. The unperturbed Hamiltonian includes 
only the on-site interaction term, which is diagonal in a num- 
ber Fock-state basis. To zeroth order the ground state |$o ) 
is the Fock state with g atoms in every lattice site. The low- 
est lying excitations correspond to the one-particle-hole (1- 
ph) states \^mn) with g + 1 particles at site m, g — 1 par- 
ticles at site 71, and exactly g particles in every other site. 
There are M{M — 1) 1-ph excitations and, because of the 
translational symmetry, they are degenerate at zeroth order 
with excitation energy U . To zeroth order the dynamic struc- 
ture factor vanishes. At first order the ground state wave- 



function is 1$ 
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J/U^2Mg{g + l)\S), where 

l^> = Ef=i(l*n«+i> + |*rm-i»/V2M is the normalized 
translationally invariant state of adjacent particle-hole exci- 
tations. In order for perturbation theory to be valid, the pa- 
rameter Jg\fM jU has to be small ]7]. This could be a sig- 
nificant restriction for systems with a large number of filled 
sites but can be perfectly realized in experiments such as 
Ref.|9j where the system has only 20 occupied sites in the 
central tube. To find first order corrections to the Af (Af — 1) 
low lying excited states we must diagonalize the kinetic en- 
ergy Hamiltonian within the 1-ph subspace. If we expand 
the eigenstates as a linear combinations of 1-ph excitations 
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En,myi«c^ml*nm) the neccssary and sufficient 



conditions that the coefficients cLv, have to fulfill are 



(5+l)(c„+lm + C„_lm)+5(c„m+l + C„TO_i) — EiC^^^, (5) 

(3) with £"■ ' = U — JEi. Besides Eq. (jSjl, the amplitudes 



have to satisfy periodic boundary conditions cj^^ 



Mm 



^nm+M 



cjj„ and the constraint c^„ 



(which prevents 
particle and hole excitations occurring at the same site). Eq. 
(|5} is analogous to the tight binding Schrodinger equation of 
a two dimensional square lattice in the xy-plane. The x direc- 
tion is associated with the position of the extra particle and the 
y direction with the position of the hole. The different weights 
g + 1 and g can be understood in the 2D-lattice model as dif- 
ferent effective masses in the two directions and the constraint 
c^„ = as a hard wall along the x — y line. The solutions 
are not straightforward due to the fact that the effective mass 
difference breaks the lattice symmetry around the x — y axis 
and makes the hard wall constraint hard to fulfill. However, in 
the limiting case of high filling factor, g ^ 1 the solutions of 
Eq. Q (including the constraints) are 
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even 



Where we used i — (r, R), with r = 1, ... M — 1 and R — 
0, . . . M — 1. The notation R' restricts the values of R to the 
ones where i? + r is an odd number and R" to the values 
where R + r is even. The constants arR — 7r(r— _R+l)/4 and 
PrR — 7r(r + i?— l + Af)/4 guarantee the orthogonality of the 
eigenmodes. 

Using Eqs. (|6} and (|8j we get an expression for the zero 
temperature dynamic structure factor given by 
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In the analytic solution, Bragg peaks have an overall enve- 
lope of the form of an inverted parabola centered at the gap 
energy U, with a maximum height proportional to sin^ {qd/2) 
and extended over an interval with an approximated width 
of 4J(2(; -|- 1) cos{qd/2). The cos{qd/2) dependence of the 
width and s\v?{qd/2) dependence of the height indicate that 
as q approaches n/d the energy transfer is highly peaked 
around U . This behavior is observed in Fig. [2 where the 
overall width decreases as q approaches Tr/d, while the peak 
height increases. In spite of the fact that the analytic solution 
uses the high filling factor approximation, the dependence on 
q of the width and peak height of the envelope is in agreement 
with the g — \ exact energy response as shown in Fig. [J 

We found no structure around hiu = 2U in the exact nu- 
merical results. This is consistent with the fact that the pa- 
rameters used in this paper lie within the regime of validity of 
first-order perturbation theory. 

For the form of energy spectroscopy we consider here, there 
is no fundamental limit to B ragg pulsed durations in contrast 
to momentum-spectroscopy llSll . However, practical consid- 
erations will likely inhibit the resolution of the individual ex- 
citation peaks in current experiments (see discussion below). 
In this case the envelope shown in Fig. ^will provide a more 
realistic depiction of the experimentally observable spectrum. 



IV. VALIDITY CONDITIONS 

The treatment we present here is based on linear response. 
In this section we indicate its strict validity conditions in terms 
of the Bragg strength Vo which is the relevant experimental 
parameter. 

After the Bragg perturbation is applied, the many-body 
state is no longer in the BHH ground state, |0). The transi- 
tion probability to an excited state \i), |ci(i)p, according to 
first order perturbation theory is given by 



Where HL^ = c":^ 



. f^R _ f,rR _ rR 



2'KqlM, and g an integer between and M — 1. The prime in 
the sum imposes the constraint g + r is even. It is important to 
emphasize that only the states with i? = have a dispersion 
relation which agrees to first order in J to the mean-field so- 
lution found in Ref. 1 10]. However, for these states iJ^ — 0. 
In Fig. [2we compare the energy transfer as a function of 
the Bragg frequency calculated from Eq. (|9} to results ob- 
tained by the exact diagonalization of the BHH for two values 
of Bragg momenta q. In contrast to the superfluid regime llal . 
where Bragg spectroscopy excites only the quasiparticle state 
with quasimomentum g, in the Mott regime we observe M — \ 
peaks pertaining to the two dimensional character of the 1 -ph 
dispersion relation. The Bragg momentum g fixes one quan- 
tum number R but the other can take M — 1 different values. 
In the analytic solution due to the constraint in Eq. (|9}, g + r 
even, only (A/— 1)/2 of the possible M—\ peaks are present. 
The constraint is a consequence of the extra symmetry intro- 
duced in the high filling factor approximation where similar 
"effective masses" are assumed. 



icwp = vimpW 



,sin^{{E,/h~uj)t/2) 
h?{E,/h-ujY '' 



(10) 



where the eigenenergies Ei of the states are measured with 
respect to the ground state energy. The validity of linear re- 
sponse requires the total excited state population at the conclu- 
sion of the Bragg perturbation to be small compared to unity: 



^ \ci{Tp 



<1. 



(11) 



Deep in the Mott regime the response of the system is domi- 
nated by the M excited states \(^\ '). Because all these states 

have energies E^ ' approximately given by C/, the maximum 
transfer energy possible is of order U . The validity of linear 
response constrains the total imparted energy to be much less 
than U and the heating rate 6E/St = U^^,^ \ci{t)\'^/t due 
to the Bragg perturbation to be much less than U/Tp. 

It was previously shown that the M excited states |<1>^ '), 
have an energy spread given by 4J(2g + l)cos{qd/2) 



> 



20 
15 
10 






(a) qd-2;r/9 






/ 


,'^ 


\ 


/ 


:;, 


\ 


/ < 


v> 


\ 


/ : 


' 


'K 


A A...^k 


■Am^* 


A jfl \ 





38 40 42 44 46 48 50 52 

120 



100 

>° 80 

::: 60 

§ 40 

20 





(c) qd-6;r/9 ^ 




/ ■ 

/ ;, 






1 ^ 

„_Ji:U*/» 




\ 



38 40 42 44 46 48 50 52 
fiu) I i 




140 
120 



38 40 42 44 46 48 50 52 
(d) qd-8;r/9 / 



^ 100 
^ 80 
W 60 
^ 40 
20 








38 40 42 44 46 48 50 52 

flCjLl I i 



FIG. 1: Energy transfer for a homogeneous system at zero temperature. Solid line exact solution, crosses perturbative solution (envelope: 
dashed line). Other parameters: M = TV = 9, .Jr/h = 20, U/J = 45. 



(1)\ 



(Eq. (|6)), and matrix elements given by: |($^ Ipjl'l'i 
Jsm{qd/2)y/32gig + l)/U, (Eq. (|9)). The average sep- 
aration between two consecutive states is of order AE ^ 
4:J{2g + l) cos{qd/2)/AI. Individually resolving the different 
lines will require one to apply the Bragg pulse for a time of 
order T^"'^ > Mh/{U{2g + 1) cos(gd/2)). The vaUdity of 
linear response, Eq. (II 1> . therefore requires that 



u fqd 

Vo«^cot - 



(12) 



For the parameters of Ref.j^l, where ^^Rb atoms are trapped 
in a lattice of depth 18.5Eji, the tunneling time h/J is about 
0.1s. The number of occupied wells is M ^ 20, with a filling 
factor 5 ~ 1. Resolving a single peak would require a Bragg 
pulse of duration Tp ' ^ 0.2s. With these conditions linear 
response is valid if Vb <C 0.015cot{qd/2)Ei^. The acceptable 
heating rate is much less than l.SEji/s. 

(s) 

If the duration of the applied perturbation is Tp < Tp , 
excited states will not be individually discernible. Near reso- 
nance (hw w U), for pulse durations smaller than the inverse 

bandwidth: Tp < tI''\ T^"^ = h/{U{2g + l) cos{qd/2)), all 



states will be resonantly excited. If Tp 
linear response requires 



14) « 



U 



cot I -— 
M \ 2 



T^^'' the vaUdity of 



(13) 



M excited states. For the parameters given above this inequal- 
ity implies Vq ^ 0.07 coi{qd/2)Eii. Here, the acceptable 
heating rate is much less than iQEn/ s. 

We note that in the superfluid regime the uncorrected na- 
ture of the system allows for a less stringent validity condition 
to hold: it is only required the amount of excited atoms to be 
small compared to the condensate population. 



V. FINITE TEMPERATURE CASE 

It is well known in the literature (see for example IVW ) that 
in the superfluid regime, Bragg spectroscopy is not an appro- 
priate tool for probing the temperature of the system. The 
reason is that even though S{q,Lo) is temperature dependent, 
experimental observables such as the energy transfer depend 
on xil, ^) which is very weakly temperature dependent. This 
is not the case deep in the Mott insulator regime. In a transla- 
tionally invariant lattice all the 1 -ph excitations have an energy 
separation of order U from the ground state and a splitting be- 
tween them of order J. If the temperature is ksT < C//3, it 
is still valid to restrict the Hilbert space to the one spanned by 
the 1-ph excitations. In this regime xili'^) (see Eq. @) can 
be written as: 



where the factor of w M, accounts for the contribution from all 
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FIG. 2: Energy transfer for a system at finite temperature. Solid fine kBT/U ■ 
same as those in FigQ 



0.007, crosses ksT/U = 0.21. The other parameters are the 



VI. FINAL REMARKS 
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(15) 

where the sum runs over the states in the one particle-hole 
band, and Z « 1 + M{M + l)e~l^^. The first term in 
Eq. (ll4t . proportional to Xg('^)|T=o, causes a thermal reduc- 
tion of the zero temperature response. The second term, which 
scales like e^^^ f]J/Z 1 19], makes the system sensitive to low 
energy Bragg perturbations at frequencies resonant with the 
energy difference between two 1-ph excitations. The factor 
e^P'^ suppresses the observability of these thermal effects for 
keT < U/b. 

In Fig. 121 we plot the energy transfer as a function of the 
Bragg frequency, as calculated from exact diagonalization of 
the BHH for two different temperatures. The figure shows that 
for temperatures kBT/U > 1/5 the height of the zero temper- 
ature peaks around U decreases, while low frequency peaks 
appear. The presence of these low-frequency peaks is there- 
fore a signature of finite temperature. In particular, if peaks 
around U are observed in absence of low frequency response, 
the temperature is lower than U/{5kB)- While this analysis 
does not provide a precise determination of the temperature, 
it is still useful because it shows that Bragg spectroscopy is 
sensitive to temperatures of order of the interaction energy. 
Current experimental techniques do not provide any informa- 
tion on the scale of U. In fact, in current experiments tempera- 
ture measurements rely on the analysis of atomic interference 
patterns after a certain time of flight following the release of 
atoms, and the measurement precision is of the order of the 
energy spacing to the second lattice band, which is typically 
one order of magnitude larger than U. 



In recent experiments [13^, Bragg spectroscopy was per- 
formed using a setup where the Bragg momentum equals the 
lattice momentum and response was observed. Our present 
analysis, in agreement with previous ones(fl(?|,fl5]), predicts 
no response for q = 2TT/d. Using similar perturbative tech- 
niques as the ones described here, we extended our calculation 
to inhomogeneous systems with a strong harmonic magnetic 
confinement. We also found no scattering for q — 2Ti/d in 
these systems. The fact that neither the inhomogeneity nor 
the finite size of the system are responsible for the observed 
signal suggests that nonlinear effects are the most plausible 
explanation for the experimental results. 

In summary, we have shown that Bragg spectroscopy can be 
a suitable experimental tool for characterizing the Mott insula- 
tor phase. By measuring the transfer energy at different Bragg 
momenta we proved it is possible to get information about 
the excitation spectrum: Bragg peaks are centered around the 
characteristic Mott excitation gap and are contained in an in- 
terval whose width is proportional to the 1-ph excitation band 
width. Their average height is maximized when the Bragg 
momentum approaches vr/d. Finally looking at the low fre- 
quency response we showed that Bragg spectroscopy is sensi- 
tive to temperatures of order of the Mott gap. 
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